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Some results about dissipative processes in superfluids are presented. 
We focus on fermionic superfluidity and restrict our analysis to the contri- 
bution of phonons to bulk viscosity, shear viscosity and thermal conduc- 
tivity. At sufficiently low temperatures phonons give the dominant contri- 
bution to the transport coefficients if all the other low energy excitation of 
the system are gapped. We first consider a system of cold fermionic atoms 
close to the unitarity limit. Then we turn to the supcrfluid phase of quark 
matter that may be realized at high baryonic density. 
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1. Introduction 

The condition for superfluidity and the hydrodynamic equations govern- 
ing its normal and superfluid components have been derived by Landau in 
his pioneering work pp. Superfluidity is a property of quantum fluids re- 
lated with the existence of low energy excitations that satisfy the Landau's 
criterion for superfluidity [H [21 [3] 

Min^l ^ o , (1) 
P 

where e(p) is the dispersion law of the excitation. 
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In general, superfluidity is due to the appearance of a condensate which 
spontaneously breaks a global symmetry of the system. As a consequence, 
in the low energy spectrum one has a Nambu-Goldstone boson, the phonon 
ip, with a linear dispersion law that satisfies Eq. (TjQ). 

Superfluidity was first discovered in 4 He, which becomes a frictionless 
fluid when cooled at temperatures below 2.17 K [21 H]. The superfluid prop- 
erty of 4 He is due to the Bose-Einstein condensation of the bosonic atoms in 
the lowest quantum state; thus quantum effects become macroscopically ob- 
servable. Fermionic systems can become superfluid as well. According with 
the Cooper theorem, fermionic superfluidity takes place in quantum degen- 
erate systems when the interaction between neutral fermions is attractive 
and the temperature is sufficiently low. In this case one has the formation 
of a di-fermion condensate that breaks a global continuous symmetry. 

We discuss the hydrodynamic equations of superfluids focusing on two 
quite interesting fermionic systems. First, we consider trapped cold atomic 
gases [5 J in the region of infinite scattering length (the so-called unitarity 
limit). Then, we turn to cold relativistic quark matter at extremely high 
baryonic densities in the color-flavor locked (CFL) phase [6j. These two 
phases of matter are quite different, however they share an important prop- 
erty: both systems are (approximately) scale invariant, and their properties 
do not depend on the detailed form of the interaction. 

The hydrodynamic equations governing the fluctuations of a superfluid 
are essentially different from standard fluid equations. In a superfluid there 
are two independent motions, one normal and the other superfluid. The 
transport properties depend on the shear viscosity coefficient, rj, on three 
independent bulk viscosity coefficients, Ci>C2>C3> an d on the thermal con- 
ductivity k. These quantities can be understood as phenomenological co- 
efficients which relate the rate of change of some quantity with the corre- 
sponding affinity [7j. The requirement that the dissipative processes lead to 
positive entropy production imposes that k, 11,(2, (3 are positive and that 
Ci — C2C3- The bulk viscosity coefficient (2 plays the role of the standard 
bulk viscosity coefficient. On the other hand, £1 and (3 provide a coupling 
between the hydrodynamic equations of the two components. The friction 
forces due to bulk viscosities can be understood as drops, with respect to 
their equilibrium values, in the main driving forces acting on the normal 
and superfluid components. These forces are given by the gradients of the 
pressure P and of the chemical potential [i. One can write in the comoving 
frame 



where P eq and (i eq are the equilibrium pressure and chemical potential, V 
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is a quantity proportional to the quantum condensate, tu^ = — (d^tp + fiu^) 
and u^ 1 is the velocity of the fluid. 

In a conformally invariant system it has been shown in Ref. [8] that 
Ci = C2 = 0. However, £3, k and r/ cannot be determined by the same 
symmetry reasoning. Regarding the shear viscosity, it is worth mentioning 
that a lower bound for the shear viscosity to entropy ratio has been derived 
employing the AdS / CFT correspondence in the strong coupling limit within 
the iV = 4 super-symmetry Yang Mills theory [9], obtaining rj/s = 1/47T. 

In the low temperature regime, TcT c , where T c is the critical tempera- 
ture for superfluidity, the transport properties of superfluids are determined 
by phonons. The contribution of other degrees of freedom is thermally sup- 
pressed. In this case one can show that Ci = C2C35 meaning that there are 
only two independent bulk viscosity coefficients and that one of the rela- 
tion for positive entropy production is saturated; the system tends toward 
the state where the velocity of the superfluid component and the velocity 
of the normal component are parallel and bulk viscosity does not lead to 
dissipation. 

For T « T c the transport coefficients strongly depend on the phonon 
dispersion law. The shear viscosity is the only transport coefficient that 
does not vanish for phonons with a linear dispersion law. But, for the bulk 
viscosities and for the thermal conductivity one has to include the term 
cubic in momentum 

e(p) = c s p + Bp 3 + 0(p 5 ). (4) 

Moreover, in the computation of the bulk viscosity one has to consider the 
processes that change the number of phonons. The parameter B determines 
whether some processes are or are not kinematically allowed. For B > 
the leading contribution comes from the Beliaev process <j) — ► (fxf). In the 
opposite case the Beliaev process is not kinematically allowed and one has 
to consider the processes cpcj) — ► cjxjxj). 

2. Cold atoms at unitarity 

Experiments with trapped cold atomic gases have reached an extremely 
high level of accuracy. The system consists of fermionic atoms, like 6 Li 
or °K, in two different hyperfine states. The fermions in the two hyper- 
fine states have opposite spin and the interaction between them can be 
tuned by means of a magnetic- field Feshbach resonance [10] . The strength 
of the interaction between atoms depends on the applied magnetic field 
and can be measured in terms of the s-wave scattering length. By varying 
the magnetic-field controlled interaction, fermionic pairing is observed to 
undergo the Bose-Einstein condensate (BEC) to Bardeen-Cooper-Schrieffer 
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(BCS) crossover. In the weak coupling BCS region the system is charac- 
terized by the formation of Cooper pairs. In the strong coupling limit the 
system can be described as a BEC dilute gas. The unitary limit is reached 
when the magnetic field is tuned at the Feshbach resonance [11], where the 
two-body scattering length diverges. 

Far from unitarity, the properties of the system are qualitatively and 
quantitatively well understood using mean field theory |12j . However, the 
mean field expansion is not reliable close to unitarity because the scattering 
length is much larger than the inter-particle distance and there is no small 
parameter in the Lagrangian to expand in. Therefore fluctuations may 
change the mean field results substantially. 

Close to the unitarity region quantitative understanding of the phases 
comes from Monte-Carlo simulations [13J, or considering the expansion in a 
small parameter that comes from the generalization to an arbitrary number 
N of spins [J3] , or considering an e = 4 — d expansion and then extrapolating 
the results to d = 3 dimensions. 

A different method comes from considering that at sufficiently low tem- 
perature the only active degrees of freedom are the phonons (see how- 
ever [IS]). The interesting point is that the effective Lagrangian for phonons 
can be determined from the pressure and by demanding non-relativistic gen- 
eral coordinate invariance and conformal invariance [16] 



C cS = com^X 5 / 2 + c.m^^t + -^(V 2 <j>) 2 Vx , (5) 

where Cq, c\ and C2 are three dimensionless and universal constants. From 
the expression above one can determine the coefficients appearing in the 
dispersion law of phonons obtaining 

and B = -TT 2 c s ^(ci + ^C2^j . (6) 
The hydrodynamic equations describing the behavior of the system are 




^ + V-(Q + Q') = 0, 

where 

Q' = q + Ki - P v n) + t ■ v n , (8) 
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and Tij, h and q are small dissipative terms that close to equilibrium are 
given by 

2 

Tij = -ri{djv ni + div nj - -SijV • v„) - <%(CiV • (p s (v s - v n )j + C2V • v„) , 

h = -C3V • (p s (v s - v n )) - CiV • v n , 
q = -kVT . 

The various transport coefficients can be determined starting from the 
phonon Lagrangian in Eq. ([5]), considering the appropriate phonon scatter- 
ing process. The shear viscosity of a unitary superfluid at low temperature 
has been computed in Ref. [T7] obtaining that 



where £ = 0.2 — 0.3 and Tp is the Fermi temperature. 

The thermal conductivity from phonons of a unitary gas has not been 
determined, yet. 

Regarding the bulk viscosity one has to consider phonon number chang- 
ing processes. For cold atoms at unitarity B > and the Beliaev process is 
kinematically allowed. In Ref. [18] it has been verified that in the conformal 
limit Ci = C2 = 0, while the remaining bulk viscosity coefficient evaluated 
within kinetic theory in the relaxation time approximation turns out to be 



The presence of non-vanishing transport coefficients leads to experimen- 
tally detectable effects. The damping of radial breathing mode depends on 
the bulk and shear viscosity coefficients. However, since Q2 vanishes, the 
bulk viscosity enters only in presence of a difference of velocity between the 
normal and the superfluid component and is in general negligible. In order 
to determine £3 one should produce oscillations where the normal and su- 
perfluid component oscillate out of phase. The transport coefficients enter 
also into the damping rate for the propagation of first and second sound 
in a superfluid [2J. The damping of first sound, a.\, depends on the shear 
viscosity and on ("2, whereas the damping of second sound, 0:2, depends on 
all the dissipative coefficients. 

3. Quark matter at extremely high baryonic density 

Relativistic superfluids might be realized in the interior of neutron stars 
where the temperature is low and the typical energy scale of particles is 



^7.7x 10- 6 £ 5 -f 
s 1 



(9) 
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extremely high. In particular, in the inner crust of neutron stars the at- 
tractive interaction between neutrons can lead to the formation of a BCS 
condensate. Moreover, if deconfined quark matter is present in the core of 
neutron stars it will very likely be in a color superconducting phase |19j . 
Quantum Chromodynamics (QCD) predicts that at asymptotically high 
densities quark matter is in the color-flavor locked phase. In this phase up, 
down and strange quarks of all three colors pair forming a difermion conden- 
sate that is antisymmetric in color and flavor indices. CFL quark matter is 
a superfluid as well, because by picking a phase its order parameter breaks 
the quark-number U(1)b symmetry spontaneously. 

There are different formulations of the non-dissipative hydrodynamical 
equations of a relativistic superfluid [201 [21] . They were derived as relativis- 
tic generalizations of Landau's two-fluid model of non-relativistic superfluid 
dynamics. The dissipative terms which enter into the relativistic hydrody- 
namical equations were derived in [22]. As it occurs in the non-relativistic 
case, for a relativistic superfluid one can define the thermal conductivity, 
k, the shear viscosity coefficient, r], and three independent bulk viscosity 
coefficients, (1^2,(3- 

The dispersion law of phonons in the CFL phase has been derived in [23] , 
and one has _ 

Co = \ — and B = ^77 . (11) 

V 3 540A 2 v ' 

In the asymptotic high density limit, the gap, A, can be computed from 

QCD [21] 

A~6 W~ 5 exp^--^ , (12) 

where b$ = 5127r 4 (|) 5//2 exp ^— ^-^j and g is the QCD gauge coupling 
constant. At very high chemical potentials, one can neglect quark masses, 
as far as m q <C [i\ moreover the coupling constant is small, g{p) <C 1, and 
one can assume that CFL is approximately scale invariant. 

A first study of the shear viscosity and of the contribution to the bulk 
viscosity coefficients due to phonons has been presented in Refs. [25] [26] . 
Beside phonons, kaons may give a sizable contribution to the transport 
coefficients. The contribution of kaons to Qi has been studied in Ref. |27j . 
There is still no computation of the contribution of kaons to the remaining 
bulk viscosity coefficients. 

Considering only the contribution of phonons one has that Ci = C2 = 
0, while the third bulk viscosity coefficient does not vanish and depends 
parametrically on the physical scales as 
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We remark that these are only approximated results that arises in the g <C 1 
limit, after neglecting the running of the QCD gauge coupling constant and 
the effect of the strange quark mass. 

The contribution to thermal conductivity due to phonons and kaons has 
recently been studied in Ref. |28j. The thermal conductivity from phonons 
turns out to be dominant and given by 



whereas the shear viscosity from phonons as determined in [25] is given by 



If superfluidity occurs in the interior of compact stars, it should be 
possible to find signatures of its presence through a variety of astrophysical 
phenomena. For example, the most natural explanation for the sudden 
spin- up of pulsars [29], the so-called glitches, relies on the existence of a 
superfluid component in the interior of the star, rotating much faster than 
the outer solid crust. After the unpinning of the superfluid vortices, there 
is a transfer of angular momentum from the interior of the star to the outer 
crust, giving rise to the the pulsar glitch. 

Another possibility to detect or discard the presence of relativistic su- 
perfluid phases consists in studying the evolution of the r-mode oscillations 
of compact stars [30j. R- modes are non-radial oscillations of the star with 
the Coriolis force acting as the restoring force. They provide a severe limita- 
tion on the rotation frequency of the star through coupling to gravitational 
radiation (GR) emission. When dissipative phenomena damp these r-modes 
the star can rotate without losing angular momentum to GR. If dissipative 
phenomena are not strong enough, these oscillations will grow exponen- 
tially and the star will keep slowing down until some dissipation mechanism 
is able to damp the r-mode oscillations. Therefore, the study of r-modes 
is useful in constraining the stellar structure and can be used to rule out 
some matter phases. For such studies it is necessary to consider in detail 
all the dissipative processes and to compute the corresponding transport 
coefficients. 




(14) 



V 




(15) 



4. Mutual friction 



Beside shear and bulk viscosity in a rotating superfluid one has to con- 
sider one more dissipative process. This is due to the scattering of phonons 
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off superfluid vortices and leads to the so-called mutual friction force be- 
tween the normal and the superfluid components. This force can be evalu- 
ated from the differential cross section per unit vortex length for the phonon- 
vortex scattering process. In case the phonons form a diluted gas, one has 
that 

da c s cos 2 9 2 

^ = M^ sm T' (16) 

where 9 is the scattering angle, E is the phonon energy and 1/A = (l — c 2 ) ^2 
with k is the quantized circulation. 



5. Addendum: Gravity analogs 

The evaluation of the cross section for the phonon- vortex scattering can 
be determined employing the gravity analogs technique. According with 
Unruh [31] , the effective action of a phonon propagating in a fluid is equiv- 
alent to the action of a scalar field in a curved background. This analogy 
can be used in two different ways; to study some aspects of general relativ- 
ity (GR), e.g. black hole evaporation, from the analogous hydrodynamical 
system; to study some properties of hydrodynamics using results of general 
relativity. In our case we want to determine the phonon-vortex cross section 
in the CFL phase, using the analogous results obtained in GR [32J. In order 
to do this we separate the Nambu-Goldstone boson field as 

<p(x) = <p(x) + (f>(x) (17) 

where <p(x) is a classical field describing the bulk motion of the superfluid 
and (j)(x) is the fluctuations on the top of the superfluid. Thus we can write 



s[ lf ] = sm + ~jd*x- 



dC, 



off 



d lx <t>d v (j> + --- (is) 



d{d^)d{d v f) 

and introducing the acoustic metric 

9iw = Viw + ( c s - 1>/^ > ( 19 ) 
we have the action for the phonon field 

S[<f>] =\( d A x^-g gTdp 4>d v <\> . (20) 



2 

When the classical background has a vortex configuration it results in a 
non-trivial metric. Then one can obtain the vortex-phonon cross section 
studying the propagation of the scalar field in this metric |33j . 
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